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ABSTRACT 

We simulate collisions involving red-giant stars in the centre of our galaxy. Such 
encounters may explain the observed paucity of highly luminous red giants within 
~ 0.2 pc of Sgr A*. The masses of the missing stars are likely to be in the range 
~ 2M — 8M . Recent models of the galactic centre cluster's density distribution 
and velocity dispersion are used to calculate two-body collision rates. In particular 
we use stellar-evolution models to calculate the number of collisions a star will have 
during different evolutionary phases. We find that the number of two-body collisions 
per star is ^ 1 in the central 0.1 pc — 0.2 pc, depending strongly on the galactocentric 
radius, with some uncertainty from the assumed cluster models and stellar-evolution 
models. Using a 3D numerical hydrodynamics code (SPH) we simulate encounters 
involving cluster stars of various masses with a 2 Mq red giant and an 8 M Q red giant. 
The instantaneous mass loss in such collisions is rarely enough to destroy either giant. 
A fraction of the collisions do, however, lead to the formation of common envelope 
systems where the impactor and giant's core are enshrouded by the envelope of the 
giant. Such systems may evolve to expel the envelope, leaving a tight binary; the 
original giant is destroyed. The fraction of collisions that produce common envelope 
systems is sensitive to the local velocity dispersion and hence galactocentric radius. 
Whereas most of our collisions lead to common envelope formation at a few parsecs 
from Sgr A*, very few collisions do so within the central 0.2 pc. Using our collision-rate 
calculations we then compute the time-scales for a giant star to suffer such a collision 
within the galactic centre. These time-scales are <; 10 9-10 years and so are longer 
than the lifetimes of stars more-massive than ~ 2 M Q . Thus the observed paucity of 
luminous giants is unlikely to be due to the formation of common envelope systems 
as a result of two-body encounters involving giant stars. 
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1 INTRODUCTION 

Recent observations of the central stellar cluster of our 
galaxy reveal that the brightest K-giant stars (2 — 8Mq 
AGB stars or more-massive supergiants) are depleted within 
~ 0.2 pc from Sgr A*, the dynamical centre of the Milky Way 
(Genzel et al. 1996). The missing stars are highly luminous 
and are therefore very large, and so Genzel et al. (1996) con- 
cluded that stellar collisions were a likely explanation for the 
paucity if the cluster's core density is sufficiently large to 
make collisions frequent and if such collisions permanently 
destroyed the photospheres of giants. 

The effects of stellar collisions on giant stars have been 
previously investigated in the context of globular clusters 
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(Davies, Benz & Hills 1991; 1992). Such collisions often re- 
sulted in the impactor becoming bound to the giant's core 
within the envelope of the giant, forming a common enve- 
lope system. In such systems, the impactor and the core 
spiral towards each other and the binding energy released 
is deposited in the envelope which is ejected on time-scales 
shorter than the evolutionary time-scale of the giant. By 
this method the giant is destroyed, leaving behind a binary 
comprised of the impactor and the giant's core. Collisions 
in galactic centres differ from those in globular clusters due 
to the higher velocities involved. The higher velocity disper- 
sions in galactic nuclei (a <; lOOkm/s) mean that incoming 
stars are faster and so are less likely to form bound systems 
than in globular cluster cores (a ~ lOkm/s). 

In this paper we investigate the frequency of two-body 
collisions involving main-sequence stars and also red giants. 
The red giants investigated are similar to those observed to 



© 0000 RAS 



2 Bailey and Davies 



be depleted in the galactic centre (hereafter GC). We simu- 
late a number of encounters involving the giants, quantifying 
the mass loss and probability of forming bound systems. If 
the GC contains binary stars, the collision time-scale be- 
tween giants and binaries may be significant if the cluster 
density is sufficiently large. In a companion paper, Davies 
et al. (1998), we extend this study to collisions involving 
binary impactors. 

In §2 we consider the stellar contents, motions and dis- 
tribution within the central parsec of our galaxy, and we dis- 
cuss the possible nature of the stars that currently cannot 
be detected. The stellar density profile and velocity disper- 
sion are then used to calculate collision time-scales in §3. 
In §4 we consider how devastating an encounter must be to 
destroy a giant. We discuss our approach to simulating col- 
lisions involving red giants in §5 and we present our results 
in §6. In a number of collisions the impactor will become 
bound to the giant star. We comment on the nature of such 
systems in §7, prior to discussing our results in §8. Alter- 
native methods of destruction for the giants are explored in 
§9. We summarise our findings in §10. 



2 STELLAR CONTENTS, DISTRIBUTION & 
DYNAMICS OF THE CENTRAL PARSEC 

2.1 Observable stars and the bright-giant paucity 

Dust in the plane of our galaxy obscures the GC by ~ 30 
magnitudes of visual extinction (Rieke, Rieke & Paul 1989). 
The infra-red K band (2.2 /im) extinction is much less severe, 
Ak ~ 3 magnitudes (Rieke, Rieke & Paul 1989). The large 
distance to the GC (~ 8.5 kpc) gives a distance modulus of 
14.6 which, coupled with Ak ~ 3, limits stellar observations 
to giant stars. Within the central ~ 0.5 pc of our galaxy 
there are ~ 600 red giants of magnitude 15 (Eckart et 
al. 1995). However, the brightest red giants are not observed 
in the deepest parts of the galactic nucleus. Observations 
show sharp drops in the numbers of 10.5 < K < 12 stars 
at ~ 0.12 pc and of K < 10.5 stars at ~ 0.2 pc (Genzel et 
al. 1996). These highly luminous stars are unlikely to have 
been overlooked by observations. Indeed, the number counts 
in Genzel et al. (1996) of 10.5 <K< 12 giants are likely to 
be 80 — 85% complete and the numbers of the K < 10.5 stars 
are 90 — 95% complete (Genzel, private communication). At 
slightly larger galactocentric radii, giants of similar luminos- 
ity to those missing have a surface-density profile, E oc p , 
for projected radius p. Assuming that this profile should ex- 
tend further into the core, we estimate the numbers of miss- 
ing giants to be twenty-five K < 10.5 magnitude stars and 
eleven 10.5 < K < 12 magnitude stars. If the surface density 
were flat within the region of missing giants then lower lim- 
its of the numbers of missing K < 10.5 and 10.5 < K < 12 
giants are eleven and six respectively. We note that the pre- 
cipitous drop, rather than a gradual one, indicates that some 
mechanism has actively destroyed at least ~ 15 — 40 stars 
in the central 0.2 pc. 

2.2 The visible stellar population 

The outcome of a collision between two stars will depend 
on their masses and sizes (i.e. evolutionary phases). It is 




Figure 1. HR diagram showing theoretical evolutionary tracks 
for stars of masses 1, 2, 3, 4, 6, 8, 10, 15, 20 Mq, models based 
upon Pols et al. (1999) with metallicity Z=0.02. The full models 
include evolution from AGB tip to remnant; we have truncated 
the tracks at the AGB tip on this diagram for clarity. The long- 
dashed lines represent contours of K of 15, 12 and 10.5, based on 
Ak = 3 and a distance modulus of 14.6 to the GC; the short- 
dashed lines are the K contours of the same values but based 
upon Ak = 4. Effective temperature is in Kelvin, luminosity is 
solar. 

useful to visualise the types of stars detectable (and those 
which are missing) at the GC using a HR diagram. For a 
star of given temperature and bolometric luminosity, its ap- 
parent K-magnitude may be calculated using bolometric cor- 
rections and V-K colours (Johnson 1966) and by assuming a 
distance modulus and a K-extinction law. Taking the galac- 
tocentric distance to be 8.5 kpc, Figure [l] shows a HR digram 
with contours of constant K appropriate for the GC, assum- 
ing Ak = 3 and Ak = 4. The effects of uncertainties in the 
GC distance on the masses and numbers of stars one expects 
to see at the GC are small compared with uncertainties in 
Ar. Rieke, Rieke & Paul (1989) report A K = 2.7; more re- 
cently Blum, Sellgren & DePoy (1996) have derived a mean 
K-extinction of 3.3 ± 0.9 but note that Ak may reach six in 
places. We see from Figure |l| that, for both values of Ak, 
the cool, visible stars must be giants and that we cannot 
easily determine their masses because of the closeness of the 
Hayashi tracks. The Figure does, however, imply that a vis- 
ible low-mass star will likely be on the AGB phase, whilst 
for stars of increasing masses the horizontal branch becomes 
increasingly visible. The duration for which a star of a given 
K-magnitude is visible is an irregular function of its mass 
(Figure ^ depending on whether either the long-lived main- 
sequence or horizontal-branch phases coincide with the cho- 
sen K-magnitude. Figure Q highlights the great sensitivity 
of the visible durations of stars on Ak- If Ak = 3, stars of 
masses <J 1.4 Mq cannot be detected at K < 10.5. If Ak = 4, 
then stars of mass 2.5 Mq cannot be detected at K < 10.5 
and the time-scale for which the stars are visible falls sharply 
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Figure 2. Visible time-scales for stars at the GC, based on stellar 
models as described in Figure The solid line denotes the time- 
scale for which the star is visible as 15 < K < 12, the short-dashed 
line represents the time-scale for which the star is visible as 12 < 
K < 10.5 and the long-dashed line shows the visible duration in 
the K < 10.5 band. Heavy lines show the time-scales for Ak = 3; 
less bold lines represent Ak = 4. The "kinks" in the profiles 
can be attributed to different parts of the stars' evolution (e.g. 
horizontal branch) overlapping with a given K-luminosity band. 
Note that stars of masses less than ~ 0.9Mq will not ascend the 
giant branch within a Hubble time. 



for masses < 4Mq. Assuming Ak = 3, stars of most masses 
are brighter than K=10.5 for about 10 6 years and most stars 
of magnitude 10.5 < K < 12 are visible for 6 - 85 x 10 5 
years. For stars of masses in the region of 6M , the hori- 
zontal branch shines in the 10.5 < K < 12 band rather than 
in the 12 < K < 15 band. We thus see a steep drop in the 
time that such stars spend in the fainter band, and a cor- 
responding peak in the brighter band. The 10.5 < K < 12 
range corresponds to the Hertzsprung Gap for stars of higher 
masses and so the 10.5 < K < 12 durations drop. The bulk 
of detectable stars sit in the 12 < K< 15 bin in which low- 
mass stars are visible for ~ 3 — 6 x 10 7 years (still as- 
suming that Ak = 3); higher- mass stars are visible for 10 7 
years. The fluctuation in visible duration for low-mass stars 
in the 12 < K< 15 range is due to changes in positioning of 
the core-helium burning phase with respect to the K = 15 
contour. Increasing Ak to four generally decreases the visi- 
ble time-scales by a factor of a few (but within an order of 
magnitude). Some stars may have their visible time-scales 
boosted in some K-magnitude range because, say, their hori- 
zontal branches now correspond to a lower K-magnitude bin 
in which, for lower Ak, only part of their giant branches had 
been detectable. 



2.3 The entire stellar population 

The contours of constant K-magnitude on Figure jl] show 
that detectable stars must be either evolved stars of a 
range of masses or very massive stars if they are on the 
main sequence. If the GC population consists mainly of low- 
mass main-sequence stars and stellar remnants then we only 
glimpse a small fraction of the total mass of the GC stellar 
population. Indeed, the detectable stars' motions and dis- 
tribution have been used to measure the total mass in the 
central parsec to be ~ 4.2 x 10 6 M , of which 2.5 x 10 6 M Q 
is the mass of the putative supermassive black hole at the 
centre of the cluster (Genzel et al. 1996, Eckart & Genzel 
1997). The amount of stellar mass in the central parsec is 
therefore about 1.7 x 10 6 M . We saw in §2.2 that only stars 
of masses ;> 1.4 Mq may be seen in the GC if Ak = 3. As- 
suming a Scalo IMF, stars of masses > 1.4 Mq can account 
for about a quarter of the mass of the population, if the mass 
contribution from stellar remnants is small. Figure ^| shows 
that low-mass giants may be detectable for about 4 x 10 7 
years, depending on the stellar mass. Given that a 1.4 M© 
star lives for ~ 4 x 10 9 years (Pols et al. 1999), for a con- 
stant star-formation rate one might expect to see ~ 3000 
stars within the central parsec. This is about a factor of two 
higher than, but is not wildly inconsistent with, the num- 
ber counts one would observe by extrapolating the observed 
numbers of stars within 0.5 pc (Eckart et al. 1995) out to a 
radius of 1 pc. 

A more detailed population synthesis was performed by 
Genzel et al. (1994); they found that the GC stellar contents 
were inconsistent with a constant star-formation history and 
Salpeter IMF. Detailed analyses are, however, hampered by 
uncertainties in Ak, possible incompleteness in the observed 
number counts and that there may not be a unique solution 
of IMF or star-formation history (or of both) which repro- 
duces the observations. 

The situation is even further complicated by the possi- 
bility of a raised low-mass cut-off to star formation at the 
GC (Morris 1993). For a cut-off of ~ 0.8 Mq and a constant 
star-formation history, one might expect roughly equal num- 
bers of white dwarfs and main-sequence stars. Higher low- 
mass cut-offs may allow stellar remnants to dominate the 
population, assuming that the population is old. A further 
complication is that the GC stellar contents may not be to- 
tally indigenous. Morris (1993) calculates the mass fractions 
of neutron stars and black holes that may have segregated 
into the central parsec from outside. Solutions computed by 
Morris (1993) allow most of the GC mass to be accounted for 
by black holes and neutron stars. These calculations, how- 
ever, suffer from requiring some parameters that are not well 
constrained, such as the low-mass cut-off and the mass of 
stellar black holes. For example, lowering the mass of stellar 
black holes from 10 Mq to 3Mq and adopting a low-mass 
cut-off of 0.08 Mq reduces the mass-fraction of neutron stars 
and stellar-mass black holes to < 10% by mass. There will 
almost certainly be non-indigenous remnants in the GC, but 
in quantities which are currently poorly constrained. 

The overall mass function of the GC population is there- 
fore very uncertain. For Ak = 3, we might expect most of 
the missing giants to be of low masses (say, 1.4 Mq — 3 Mq), 
for most IMFs. However, the uncertainty in Ak could mean 
that the brightest missing stars are of masses <; 4 M . 
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2.4 Distribution and motions of GC stars 

An observed number density of stars traces that of the whole 
cluster if the mean mass of the cluster stars is close to the 
mean mass of the observed stars and if the system is relaxed. 
The relaxation time (see Binney & Tremaine 1987) at 0.2 pc 
is ~ 2 x 10 9 years for a 1 Mq star, so we expect low-mass 
(1 — 2 Mo) giants to trace the overall distribution of stars 
of similar masses. Relaxation time decreases with the mass 
of the star but more-massive stars do not live long enough 
to become relaxed in the GC. Such stars should be treated 
cautiously as mass tracers. Genzel et al. (1996) model the 
observed red-giant number-density distribution to be 



n(r) oc 



l + (r/4.2") : 



(1) 



for galactocentric radius, r. We adopt this Equation, as- 
suming that the observed stars are indeed tracing the whole 
population. This assumption is somewhat supported by the 
stars' isotropic motions (Genzel et al. 1996, Eckart & Genzel 
1997), indicating a relaxed population. Genzel et al. (1996) 
also give the central stellar density, po, to be 4 x 10 6 Mq/pc 3 
and a core radius of 0.38 pc. We normalise Equation (^) with 
this value of po but note that po and the core radius are inter- 
linked values and that the latter is subject to some debate, 
depending on whether the cluster's surface-brightness distri- 
bution or the surface-density distribution is used (see Gen- 
zel et al. 1994). By varying po and the core radius between 
0.25 pc and 0.5 pc whilst conserving the total central-parsec 
mass within errors set by Genzel et al. (1996), we find that 
the collision time-scales for stars (see §3) can vary by up to 
a factor of 1.5. 

We also note that Equation ([[]) gives quite a flat pro- 
file inside the core radius. However, observations become 
crowded within this radius and it becomes harder to con- 
strain the density distribution. If the condensed mass at Sgr 
A* is a black hole, theoretical modelling predicts a density 
cusp at low radii: p oc r ' where 7 ~ 1.5 — 2.5 (Quin- 
lan, Hernquist & Sigurdsson 1995; see also Bahcall & Wolfe 
1976). Consequently we also calculate collision rates assum- 
ing p oc r~ 1,8 (coinciding with the higher-radius density pro- 
file at the GC), which we set equal to the density profile of 
Genzel et al. (1996) for r > 0.38 pc where we expect the 
latter to be accurate. 

A final comment on Equation (jl]) is that it is based 
upon (M/L)k = 2, the mass-to-light ratio evaluated at a 
few parsecs from the centre. Genzel et al. (1996) claim that 
their usage of only brighter cluster members reduces the 
effects of population change at lower radii. There does exist 
some evidence of an increase in (M/L)k at low radii: Saha 
et al. (1996) show that this ratio may increase by a factor 
of 3-5 between 2.5 pc and 0.5 pc, although if one subtracts 
2.5 x 10 6 Mq for the putative central black hole then the 
increase in (M/L)k for the cluster only reduces to two. 

Collision-rate calculations require a velocity dispersion 
(see §3). We adopt the deprojected velocity dispersion given 
by Genzel et al. (1996): 



55 2 + 103 2 (r/10") -1 ' 2 (km/s) 2 



(2) 



The velocity dispersion rises as galactocentric radius falls 
due to the presence of the supermassive black hole (see also 
Eckart & Genzel 1997) . At higher radii, the enclosed cluster 



mass dominates and a flattens out. Equation (Q) is valid 
throughout the central parsec. 



3 COLLISION RATES OF THE 
CENTRAL-PARSEC STARS 

The collision rate per star of species 1 with stars of species 
2 is 

ri = £i2(r) = f £ £y Rmix^VaoP^Voa) dVoo (3) 

"iW Jo 

where n is the number density of stars and P(Voa) is the 
probability that two stars have relative speed Voo ■ S g is the 
gravitationally focused cross section: 



E g = 7ri? min 1 + 



2 (, , 2 G(Mi + M 2 ) 

R min Voo 



(4) 



where Mi and M2 are the stellar masses and i? m i n is the 
collision radius (which we take to be equal to the sum of the 
radii of the two colliding stars) . The second term on the right 
in Eg represents the mutual attraction between two stars 
which increases the cross section above that of the geomet- 
ric value. We take P(Voo) to be Maxwellian, which is appro- 
priate at high galactocentric radii where the stellar cluster 
dominates the enclosed mass. In the region where the black 
hole dominates the enclosed mass (say, <J 0.7 pc, see Gen- 
zel et al. 1996, Eckart & Genzel 1997), P(Voo) is still close 
to Maxwellian (Quinlan et al. 1995). At any galactocentric 
radius we take the relative velocity dispersion between two 
stellar species of any mass to be 2cr 2 , where the variation of 
a with galactocentric radius is given by Equation (|). In a 
system of stars that has achieved equipartition of energy, we 
have a 2 oc 1/M (see Binney & Tremaine 1987); stars of dif- 
ferent masses would have different velocity dispersions. The 
effects of equipartition on the collision rates do not change 
our conclusions. 

Assuming number- density and velocity-dispersion pro- 
files discussed in §2, Figure ^ shows collision time-scales for 
the GC. The collision time-scale is a strong function of r and 
is shortest deep within the core. Comparing Figure ^| with 
Figure |^ we see that the time-scales for collisions involving 
giants are similar to the durations for which they are visible 
at the GC. The collision time-scale for solar- type stars is 
considerably longer. It is only noticeably less than a Hubble 
time at very low radii if the density profile is a cusp (i.e. 
within ~ 0.1 pc, assuming p oc r~ 1,8 ). 

The calculations shown in Figure [| have used an in- 
stantaneous radius for the stars. Such calculations lose their 
meaning, however, if the stellar radius changes on time- 
scales shorter than the collision time-scale. We now calculate 
the number of collisions that a star will have by integrating 
the instantaneous collision rate over its lifetime: 



n-coii = Vidt 



(5) 



In this way we obtain Table [l] in which we have re- 
stricted r to interesting parts of the star's evolution as well 
as the entire lifetime, for stars of various masses. We see 
that n C oii, when evaluated over the lifetime of a star, can be 
a substantial fraction of unity. For sufficiently low radii, and 
preferably with a steep cusp density profile, n co n is about 
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Mass [M e ] r [pc] 


Profile 










n coll 












LIFE 




AIS 




K< 15 (A K = 3) 


K< 15 (A K = 4) 


1 0.12 


G96 


0.741 




0.292 




0.278 




0.229 




0.12 


cusp 


1.74 




0.686 












0.17 


G96 


0.549 




0.242 




nisi 

U. ±ol 




n i As 




0.17 


cusp 


0.894 




0.394 








n 949 




1.5 0.12 


G96 


0.550 




0.129 




0.320 




0.272 




0.12 


cusp 


1.29 




0.302 




U. ( 01 




U.D'rtU 




0.17 


G96 


0.387 




0.106 




n 9ns 
u.zuo 




n 1 7fi 




0.17 


cusp 


0.631 




0.172 




u.ooy 








2 0.12 


G96 


0.549 




0.077 




0.340 




0.313 




0.12 


cusp 


1.29 




0.182 




n 7Q7 




n 7^a 




0.17 


G96 


0.379 




0.0636 




n 91 q 




n 90i 




0.17 


cusp 


0.619 




0.104 




U.OOo 




n ^97 




4 0.12 


G96 


0.533 




0.0246 




0.510 




0.508 




0.12 


cusp 


1.25 




0.0576 




1.19 




1.19 




0.17 


G96 


0.354 




0.0205 




0.333 




0.333 




0.17 


cusp 


0.576 




0.0335 




0.524 




0.542 




8 0.12 


G96 


0.909 




0.0116 




0.900 




0.898 




0.12 


cusp 


2.13 




0.0272 




2.11 




2.11 




0.17 


G96 


0.589 




0.00990 


0.582 




0.580 




0.17 


cusp 


0.961 




0.0161 




0.948 




0.945 




Durations [yr] 


1 




1.25 x 


10 10 


1.10 x 


10 10 


5.22 x 


10 7 


1.78 x 


10 7 


1.5 




3.10 x 


10 9 


2.73 x 


10 9 


4.69 x 


10 7 


1.09 x 


10 7 


2 




1.50 x 


10 9 


1.17 x 


10 9 


2.73 x 


10 7 


4.53 x 


10 6 


4 




2.15 x 


10 x 


1.79 x 


10 8 


3.56 x 


10 7 


3.52 x 


10 7 


8 




4.25 x 


10 7 


3.72 x 


10 7 


9.61 x 


10 6 


5.21 x 


10 6 



Table 1. Phase durations and number of collisions, n co n, during various evolutionary phases for stars of masses IMq, 1.5 Mq, 2Mq, 
4Mq and 8Mq at galactocentric radii of 0.12pc and 0.17pc. Profile denotes stellar number-density profile used: G96 refers the density 
of Genzel et al. (1996) and cusp refers to a density profile of p oc r~ , set equal to its counterpart at r > 0.38 pc. We have assumed that 
the mean mass of the GC stars, M2, is IMq but note that, approximately, n co vj oc I/M2. LIFE represents the life-span of the star from 
ZAMS to AGB termination, MS represents the main-sequence duration and K < 15 represents the time for which the stars are brighter 
than K=15 for a given extinction. Evolutionary models of stars are from Z=0.02 models from Pols et al. (1999). 

mass-loss rate by a factor of 10 in the models of Pols et al. 
(1999) produces a factor of two change in n co n. 

The relative likelihood of collisions during different 
phases of the stars' evolution are shown in Figure ^| We 
see that one of the stars in single-single stellar collisions 
in the GC will most likely be a giant star (provided that 
it is sufficiently old to have turned off the main sequence). 
In terms of collision frequency, the AGB becomes increas- 
ingly predominant over the first giant branch for stars of 
higher mass. Only a minority fraction of target stars are 
on the main sequence when hit by an impactor. Although 
the main-sequence phase is much longer lived than the giant 
phase, more two-body collisions in the GC will occur whilst 
the target star is a giant because of its much greater radius 



one. The more-massive stars are more luminous than the 
1 Mq star and so spend a larger fraction of their lives as 
K < 15 magnitude stars. During this period they reach larger 
radii than 1 Mq stars (Figure hi) and so we expect them to 
undergo greater proportions of collisions as detectable gi- 
ants. There are a number of uncertainties in n C oll- Uncer- 
tainties in the stellar number-density profile are discussed 
in §2.4. The value of the mean stellar mass (M2) is uncer- 
tain. In Table |l| we set the number of impacting stars to be 
Ti2(r) = pon(r)/M2, where n(r) is given by Equation (|l|), 
and where M 2 = 1 M . Thus, via Equation ([|, Ti oc I/M3 
(ignoring gravitational focusing). A decrease in M2 to, say, 
0.7Mq, increases Fi by a factor ~ 1.4. n co li is also sensi- 
tive to the radius-time behaviour of the stars. Changing the 
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Figure 3. Collision time-scales (in years) for stars in the GC as a 
function of galactocentric radius (in parsecs). For the calculation 
we have assumed that the background cluster members are of 
mass 1 Mq , Equation (0) for the velocity dispersion, and for the 
density profile we have used Equation (1) (solid lines) or p oc r~ 1,8 
(dashed lines), set equal to Equation (I) for r > 0.38 pc. The two 
lower lines represent the time-scale for encounters for 8 Mq stars 
at R m i n = 243 Rq . The middle two curves represent the time- 
scale between collisions for 2 Mq stars at i? m i n = 94 Rq . For 
comparison, the upper-most lines show the encounter time-scale 
profile for a 1 Mq star at i? m i n = 2 Rq. 



Figure 4. The distribution of collisions, n co n (normalised), dur- 
ing the lifetimes of stars of various masses, plotted as a function 
of the stellar radius (solar units), which emphasises different evo- 
lutionary stages. The solid line represents a 1 Mq star, the short- 
dashed curve represents a 1.5 Mq star, the long-dashed curve cor- 
responds to a 2Mq star, the short-dashed, dotted line represents 
a 4Mq star and the long-dashed, dotted curve corresponds to 
an 8Mq star. We have assumed that all impinging stars are so- 
lar mass and solar radius, and have adopted a = 183.13 km/s, 
appropriate for a galactocentric radius of 0.17 pc. 



compared to the main- sequence radius. In fact the number 
of target stars on the main sequence is <J 25%, with the ex- 
ception of 45% for a 1 Mq star. The importance of the main- 
sequence phase decreases with increasing stellar mass. The 
horizontal branch also generally dominates over the main 
sequence in competition for collisions: it is shorter-lived but 
the star is considerably larger on the horizontal branch than 
when on the main sequence. 



4 CRITERIA FOR DESTRUCTION OF 
GIANTS 

To determine whether stellar collisions may destroy giant 
stars, one needs to know how much envelope mass a gi- 
ant must lose before it is destroyed. Strictly speaking, to 
satisfy the observed paucity we must strongly reduce the 
K-magnitude of the brightest giants. Considering mild en- 
counters where the mass lost is insufficient to alter the gi- 
ant's long-term evolution, we actually expect a transient in- 
crease in the stellar K-luminosity: the impactor loses some 
of its kinetic energy to the envelope which consequently ex- 
pands somewhat, becomes redder and therefore brighter in 
the K-band. The envelope will then settle back to normal 
on its thermal time-scale (~ few- hundred years). To explain 
the paucity of giants, collisions must therefore destroy the 
giant, i.e. they must cause large amounts of envelope loss. 
The actual amount of mass loss required as a direct result 



of a collision is uncertain but is likely to be most if not 
nearly all of the envelope (Eggleton, Tout, private commu- 
nications). A second way to destroy the giant is via indirect 
mass loss after a collision produces a common envelope sys- 
tem. In such systems, the binding energy released by the 
impactor and core spiralling in towards each other may un- 
bind the envelope. To form a common envelope system in a 
hyperbolic collision requires the impactor to deposit at least 
[iV^o /2 of its orbital energy into the giant but need not nec- 
essarily produce copious mass loss. Thus to destroy a giant, 
a collision must involve either immediate, large mass losses 
or the formation of a common envelope system which then 
evolves to expel the envelope (§7). 



5 APPROACH TO HYDRO DYNAMICAL 
SIMULATIONS 

We simulated encounters involving giant stars using a 
smoothed-particle hydrodynamics (SPH) code. For a discus- 
sion of SPH, see Benz (1990). The impactors encountering 
the star were point masses which were adequate to simu- 
late either main-sequence stars or compact remnants such 
as white dwarfs, neutron stars or stellar black holes because 
of the size contrast between these objects and the giant's 
envelope. 
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5.1 Collisions performed 

We performed collisions on two red giants. The radii and lu- 
minosity of the stars were selected so that their positions on 
a HR diagram coincided with those of the missing bright gi- 
ants (see Figure [l]). One giant had mass 2Mq, radius 94 Rq 
and core mass 0.52 Mq. The other giant had mass 8Mq, 
radius 243 R0 and core mass 2.5 Mq. In each case we repre- 
sented the core by a point mass. The masses of the two stars 
were selected so that we sampled both ends of the likely mass 
range of the missing giants (see §2). Both giants were based 
upon models of temperature and density profiles discussed 
in Pols et at. (1995). 

Most simulations were performed with ~ 8500 parti- 
cles. A number of simulations were run with 16 x 10 3 and 
36 x 10 3 particles in order to quantify resolution effects 
(see §6.3). A wide range of impactor masses were selected. 
For the 2 Mq giant we used 0.6 M©, 1M©, 1.4M©, 2M , 
4 M© and 8 M© impactor masses. Our results suggested that 
for the 8M© giant, results from impacts with a 1M© im- 
pactor would be quantitatively similar to the outcomes of 
the giant's encounters with 0.6 M© and 1.4 M© impactors, 
so for the more-massive giant we simulated collisions involv- 
ing only 1M©, 2M©, 4M© and 8M© impactors. 

For each giant-impactor pair, we ran ~ 16 — 25 simu- 
lations distributed on a Voo — -Rmin grid. A number of simu- 
lations were also performed very near the Voo and i? m i n axes 
in order to quantify our extrapolations from the grid to low 
Voo and low R m i n values (see §6.3). In total approximately 
200 runs were performed. 



500 



-500 




500 



-500 

-1 I 1 1 1 1 I 1 1 1 

500 -500 
x 

Figure 5. Time sequence of fluid-particle plots of an SPH simula- 



tion between an i 
frame), with Voo 
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giant and a 1 M© impactor (centre-of-mass 



50km/s and R n 



100 R©. Here we take a 



slice in the plane of the collision. The trajectory of the impactor 
is clearly deflected; this is due to the giant's mass enclosed in its 
orbit. The top-left and top-right portions in particular show the 
shock caused by the impact. Units of x and y are solar radii. 



5.2 Quantifying collision outcomes 

To determine whether a collision destroyed the giant, we 
measured the mass loss, AM, searching for mass losses of 
most of the envelope (§4). To determine AM, the enthalpy 
of each particle was evaluated relative to both the impactor 
and the giant, using a similar method as prescribed in Ra- 
sio & Shapiro (1991). For each collision we ascertained the 
amount of mass remaining bound to the giant, the amount 
that had become bound to the impactor and the amount of 
escaping matter. We found that AM settled down within a 
sound-crossing time or two after a collision. 

A number of collisions involved dissipation of enough 
of the impactor's orbital energy for the impactor to be- 
come bound to the giant. Such systems are important as 
they potentially lead to the destruction of the giant via a 
common envelope phase (§6.2). To delimit the region of the 
Voo — Rmin plane in which two impactors became bound, 
we measure the fractional change in their orbital energy, 
AE/E . We define AE such that AE / E = 1 when the im- 
pactor just becomes bound to the giant. The initial orbital 
energy is simply the final orbital energy is calcu- 

lated from the masses, positions and velocities of the two 
centres of mass. 

For some of the bound systems formed, the impactor re- 
mained inside the envelope, i.e. a common envelope system 
had effectively been formed immediately. In such systems 
we do not have two distinct objects to which we can assign 
fluid particles. Our algorithm is thus unsuitable for such out- 
comes and so AM and AE/E were not measured for these 
outcomes. 



6 RESULTS 

The impactor hits the giant supersonically, shocking the en- 
velope, as can be seen in Figure |B| Due to the giant's grav- 
ity, the impactor is pulled around the core somewhat; the 
amount of deflection depending on the Voo of the encounter. 
Almost all of the escaping particles had spee ds in excess 
of the local escape speed (Vesc = ^JlGM^/R for R > Ri, 
where Ri is the giant's radius), demonstrating that the es- 
caping fluid was energetically unbound by virtue of its ki- 
netic energy rather than its thermal energy. Some fluid had 
speeds up to four times that required to escape and thus car- 
ried off excess energy thereby reducing the mass loss. Given 
that mass losses are generally low (§6.1), this implies that 
deep collisions are required to remove envelope material lo- 
cated close to the core. 

Both mass loss from the giant and orbital energy loss 
from the two stars increase with increasing impactor mass 
and with decreasing Voo , indicating the dominance of time- 
dependent gravity over shocks as the energy-transfer mech- 
anism. Mass loss and orbital energy loss also increase with 
decreasing 7? m i n . 

We divide the collision outcomes into two categories: 
bound and unbound giant-impactor pairs (after collision). 
A bound system is defined by AE/E > 1 (see §5.2). Both 
mass loss (Figure ^) and orbital energy loss seemed to vary 
smoothly across the boundary (Figure . 

6.1 Mass loss 

We reiterate that the instantaneous mass loss increases with 
impactor mass and falls with increasing Voo and -Rmin- For 
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two given colliding species, this implies that we would ob- 
serve maximal mass loss for collisions of lowest Vx> and low- 
est Rmin. Those particular encounters, however, tended to 
directly produce common envelope systems (see §6.2) for 
which we do not measure the mass loss due to the collision. 
Consequently the subsequent comments concerning AM are 
limited to outcomes that are not common envelope systems 
and that constitute the larger portion of our data. 

Collisions of both giants with impactors of mass ~ 
IMq involved small mass losses, <J 0.1 Mq (and usually 
much less, see Figure []). In fact, for the 8Mq giant and 
for Voo <; lOOkm/s, the mass losses were so small that too 
few SPH particles escaped to take a reliable measurement of 
AM. Even deep impacts (R m i n ~ Ri /4) were not devastat- 
ing for either giant, although the mass lost from the 2 Mq 
giant rose to ~ 0.1 Mq. 

If the bulk of the population in the GC are stars of 
masses <^1Mq, then the paucity of bright giants cannot be 
explained by instantaneous mass loss due to two-body colli- 
sions because the vast majority of collisions with 1 Mq im- 
pactors will involve mass losses much less than the envelope 
mass. 

Mass losses were higher for more-massive impactors be- 
cause of the greater gravitational forces involved. We dis- 
cuss here the effects of 8Mq impactors. Encounters involv- 
ing our 8Mq giant showed mass losses of ~ 0.1 — IMq 
for i? m in <; 50 Rq. Collisions involving the 2Mq giant with 
i?min > 25 Rq gave mass losses ^ 0.5 Mq, dropping to a 
few x10~ 2 Mq for fast (250 — 450km/s) grazing encoun- 
ters. Deeper impacts had much stronger effects for both gi- 
ants. For our 8Mq giant, collisions with R m in < 50 Rq gave 
rise to mass losses of about 3Mq for Vx> ^ 150km/s (our 
slower collisions gave rise to common envelope systems) . For 
the 2 Mo giant, encounters of R m in ,Sa 25 Rq involved mass 
losses approaching the entire envelope mass. It is possible 
that low i?min (is 30 Rq) collisions with massive impactors 
could destroy the giants. Such collisions, however, will be 
rare; considering the low -Rmin alone places upper limits of 
frequencies of 5% of all giant-8 Mq impactor collisions for 
the 8Mq giant and 20% for the 2Mq giant, assuming that 
all collisions occur at Voo = V2a(r) and that r — 0.2 pc. 
We note that in §3 we showed that a star could expect at 
most about one collision before becoming a remnant; most 
collisions will occur at too high -Rmin to destroy the giant. 
Moreover, to destroy all giants in this manner would further 
require the GC to contain a large proportion of impactors 
of masses <; 8 Mq . Such a population seems quite unlikely, 
unless, perhaps, the GC core is dominated by massive black 
holes which have diffused inwards from higher galactocentric 
radii (see §2 and Morris 1993). 

In summary, mass losses due to most two-body collisions 
are likely to be small (sJ 0.1 Mq), and though very massive 
impactors might be able to destroy the giants, the deep- 
envelope collisions required will make such encounters rare, 
unless the GC contains an exotic population of numerous 
massive stars which is itself unlikely. 

6.2 Formation of bound systems 

For collisions of sufficiently low Voo and Rmin, the impactor 
lost enough kinetic energy to become bound to the giant. 
Such collisions will almost always be physical collisions and 
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Figure 6. Impactor-accreted mass (unbroken circles) and mass 
lost from the system in collisions (broken circles) between the 
2 Mq giant and a 1 Mq impactor, as a function of Voo [km/s] and 
^min [Rq] - The symbol "B" represents that the collision product 
is a binary, and "CE SYSTEM" denotes the direct formation of a 
common envelope system. The size of the circles is proportional 
to the log of the mass. Masses are in solar units. 



we expect the resulting orbits' periastrons to remain inside 
the giant. If the apastron of the orbits were also within 
the giant's radius then a common envelope system had ef- 
fectively immediately been formed. Systems with apastrons 
outside the giant's envelope were very eccentric and had pe- 
riods of a few years to centuries. Thus the time-scale for 
further collisions is much shorter than the giant's evolu- 
tion time-scale. We expect each subsequent encounter to 
dissipate more of the orbital energy until the apastron falls 
within the giant's radius (or the envelope puffs up enough to 
absorb the impactor's orbit), producing a common envelope 
system. 

Common envelope systems are of great interest here 
as their evolution may quickly destroy the giant. The im- 
pactor and giant's core spiral towards each other, transfer- 
ring angular momentum and energy to the envelope. If a 
sufficient quantity of the potential energy between the im- 
pactor and giant's core is deposited into the envelope, the 
entire envelope may be expelled, leaving a tight binary. Al- 
ternatively, the impactor may merge with the core before 
envelope ejection is complete, terminating the common en- 
velope phase. The outcomes of common envelope evolution 
involving either of the two giants considered here are dis- 
cussed in §7, here we calculate the fraction of collisions that 
produce them. 

Figure ^ shows the dependence of the fractional change 
in pre-collision orbital energy on Vx> and R m i n for a cer- 
tain giant-impactor pair. Bound systems are produced in 
the bottom-left of the Figure and are delimited from un- 
bound systems by the AE/E — 1 contour (which we now 
denote as Rb(Voo))- As impactor mass increases, the de- 
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Figure 7. AE/E contour plot for collisions involving the 2Mq 
giant and a 1 Mq impactor, as a function of Voo [km/s] and R m i n 
[Rq]. The contours are logarithmically spaced by 0.5, the bold 
line represents AE/E = 1 i.e. the threshold of bound-system 
formation. Contours to the bottom- left of this line have AE/E > 
1 so the bottom-left region represents the region where bound 
systems are formed. Collisions of the upper-right portion leave 
the two stars unbound. 



limiting contour moves toward the top-right portion of the 
Figure. Values of AE/E between data points are obtained 
by linear interpolation of a fitted 2D exponential function 
AE / E(Voo, -Rmm). Note we do not measure AE/E for com- 
mon envelope systems, rather, we extrapolate from the other 
data. The fraction of collisions, / ce , that fall underneath the 
■Rb(Voo) curve depends on the likelihood of having a given 
velocity (Voo) and therefore on P(Voo), the velocity distri- 
bution function, which we take to be Maxwellian (§3). For 
a given Voo , fee depends on the proportion of collisions that 
have i? m in < -Rb(Voo) and thus depends on the amount of 
gravitational focusing (Equation Q): 

_ E g (Voo , flmin = i?b ( Voo ) ) VooP( Voo ) dVx, 

_ £g(Voo,iCin = i?l)VooP(Voo)dVoo 

The dependence of a(r) on position within the GC exerts 
a galactocentric radial dependence on / cc and is shown in 
Figure |§| 

We see that / co increases with impactor mass, which 
follows because the increased gravitational interaction al- 
lows bound systems to be formed at both higher speeds 
and at higher i? m i n . Figure ^ shows that for increasing V», 
the largest i? m i n that produces a bound system decreases. 
As galactocentric radius decreases, the median Vx> increases 
and so we expect / CG to decrease rapidly inside the core ra- 
dius (~ 0.38 pc) and to flatten at high radii in accordance 
with the nature of the velocity dispersion (Equation as 
illustrated by Figure ^. Within the core radius the fraction 
of bound systems formed is only a few percent (~ 1 — 10% 
at the core radius). At a few parsecs, / ce ~ 10 — 40% for 



Figure 8. The fraction of encounters that produce bound systems 
as a function of galactocentric radius (in parsecs): the 8Mq giant 
(lower half) and the 2Mq giant (upper half). Different curves 
reflect different impactor masses with curves of increasing height 
representing increasing impactor mass. For the 8Mq giant, the 
impactor masses are IMq, 2Mq, 4Mq and 8Mq. The same 
notation has been used for the 2Mq giant, except we also show 
curves for 0.6 Mq impactors (dotted line) and 1.4 Mq impactors 
(dot-dashed line). We have used Equation (^) to describe <r(r). 

the range of impactor masses considered. The large / ce at 
high galactocentric radii is due to the low local velocity dis- 
persion. 

We note that / co is substantially less than one, even for 
collisions between the 8Mq giant with 8Mq impactors at 
high galactocentric radii (i.e. low speeds). Since / ce is the 
fraction of physical collisions that lead to the capture of the 
impactor, for tidal capture to occur we require / cc > 1. Thus 
tidal capture by giants similar to those discussed here will 
be extremely rare in the galactic centre. 

6.3 Resolution and error considerations 

Three sources of error in our data were considered: the ac- 
curacy of extrapolation of AE/E to regions near the Voo 
and i?mm axes (Figure ^), errors in energy conservation in 
individual SPH runs and errors due to resolution effects (i.e. 
number of particles used) in our simulations. 

To assess errors in AE/E extrapolation we performed 
a number of simulations near the Voo and i? m i n axes and 
compared the local extrapolated AE/E value to that ob- 
tained from the extra runs. We then locally refitted -Ri(Voo) 
and re-calculated / cc to evaluate its error. Errors in / ce due 
to energy conservation in individual runs were examined by 
estimating their effects on the local AE/E value and hence 
on the locus of the -Rb(Voo) contour, producing an error in 
/ ce . In both cases the error in / co for all giant-impactor pairs 
throughout most of the GC is within ±30%. Slightly over- 
widely spaced grids for the 2 Mq giant and 2 Mq , 4 Mq and 
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8Mq impactors led to larger errors at low galactocentric 
radii, but the net effect is likely to be a mild overestimation 
(within a factor of two) of / ce for any galactocentric radius 
outside of the central few hundredths of a parsec. These er- 
rors do not weaken our conclusions; indeed, overestimation 
of fee actually acts to strengthen our conclusions (see §8). 

We also estimated errors in AA1 and AE/E (and hence 
/ ce ) due to resolution, i.e. using a limited number of parti- 
cles in our simulations. We repeated a quantity of collisions 
involving the 8M giant with 16 x 10 3 and 36 x 10 3 par- 
ticles. The quantity of interest (e.g. AM) was evaluated as 
a function of the particle separation. Extrapolation to zero 
spacing (i.e. infinite resolution) obtains the best estimate of 
the quantity. We found that our mass losses became gen- 
erally good estimates for AM <; few x 0.01 — 0.1 Mq, de- 
pending on impactor mass. The resolution errors in AE/E 
caused uncertainties in / ce that were generally small and 
never large enough to affect our conclusions. Since the 2 M 
giant shows, to some extent, similar scaling to the 8 M 
giant (e.g. number of particles per fractional radius), we ex- 
pect similar resolution errors. 

The uncertainties considered here do not affect our con- 
clusions, and so we have not included them in our calcula- 
tions. 



7 EVOLUTION OF BINARY SYSTEMS 
FORMED 

We now consider the outcomes of the evolution of com- 
mon envelope systems formed by our two-body encounters. 
Within the envelope, the core and impactor orbit each other 
in ever tightening orbits, their potential energy being con- 
verted into their orbital kinetic energy and thence deposited 
in the envelope. If the envelope gains enough energy, it will 
become unbound. We take the binding energy of the enve- 
lope to be (de Kool 1990): 

G(M C + Mmv)M e 



(7) 



where M c is the giant's core mass, M env and R env are respec- 
tively the envelope's mass and radius and we take A to be 
0.5. The change in binding energy between the impactor and 
the giant from the onset of common envelope evolution to 
its termination when the entire envelope has been expelled 
is given by: 



GM 2 M C GM 2 (M C + M m 



(8) 



2d f 2di 

where M2 is the impactor mass, di and d{ axe the initial 
and final semi-major axes. By assuming all the envelope is 
removed with some efficiency, a, then 



E c 



aAE ti . 



(9) 



Thus we may calculate the final separation between the im- 
pactor and core: 



di = M 



/ 2E C 



+ ■ 



M c + M c ,_ 



(10) 



V GaM 2 ' di 

The following results are insensitive to di since di must be 
very deep in the envelope before the second right-hand term 
becomes comparable to the first. Numerical modelling of 
common envelope evolution suggests that a ~ 0.24 — 0.6 



(Taam & Bodenheimer 1989, Sandquist et al. 1998); we take 
a — 0.4. Using the parameters for our giants (see §5), for 
both giants we find that df takes values between ~ 0.8—4 R 
for impactor masses in the range 0.6 Mq to 3Mq. df is sim- 
ilar for both giants because the greater energy required to 
unbind the 8 M giant's envelope is offset by the greater en- 
ergy released by an impactor falling towards its more mas- 
sive core. 

We have thus far assumed that the core and impactor 
are point masses. Their sizes, however, can be important 
and can lead to termination of the common envelope phase 
before total ejection of the envelope. Either the impactor 
or the core may fill their Roche lobes or, alternatively, the 
impactor may fall so deep into the giant that the local den- 
sity is comparable to its own. In the following calculations 
we find the former almost always occurs before the latter. 
Either way, the impactor is no longer distinct from its sur- 
roundings and no further potential energy may be liberated. 
The occurrence of such "premature" termination depends on 
the size of the impactor and the other parameters in Equa- 
tions 0-0. We ascertain whether envelope ejection will occur 
in the following way. Two species of impactor are considered: 
ZAMS stars and stellar remnants. The radii of the former 
are taken from Tout et al. (1996). A simple mass-radius rela- 

— 1/3 

tion for white dwarfs suffices for this discussion: R oc M x 
(Mi < 1.4 M ), scaled to 0.0248 R at 1M . More-massive 
remnants (neutron stars) are given a size of 10km. We cal- 
culate the impactor's Roche-lobe size using Eggleton (1983): 



(11) 



Rli _ 0.49g3 

d f 0.6qi + ln(l + qi) 

where q = MijM c . Finally, we used the density profiles of 
our giant models (see §5) to determine the giant's radius at 
which the density was comparable to that of the impactor's 
density (derived from its mass and the appropriate mass- 
radius relation) . Our results are shown in Figure ^[ 

We see that for a ~ 0.5, ZAMS stars of masses <; 1 M 
may expel the envelope of the 2M giant, whereas ZAMS 
stellar masses in excess of 2M are required to expel the 
8M giant's envelope. Increasing a decreases the minimum 
mass required. We also see that a ZAMS star will always fill 
its Roche lobe before penetrating sufficiently deep into the 
envelope for the local density to be comparable to its own. 
Stellar remnants should almost always be able to eject the 
envelopes of either star, the exception being for low-mass 
(^ 0.4 M ) remnants with the 8M giant. If a is very low, 
the core will generally fill its Roche lobe before a remnant 
counterpart does, except in the case of the 2 M giant with 
remnants of masses <J 0.5 M . This is because the 2M 
giant's core (of mass 0.52 M ) is essentially a white dwarf 
and white dwarfs of lower masses are larger: they will fill 
their Roche lobes first. A similar effect is not seen for the 
8M giant: its core is much larger (0.25 R ) than a stellar 
remnant and tends to fill its Roche lobe first. 

It seems likely that most low-mass main-sequence stars 
will not expel either of the giants' envelopes, particularly 
that of the 8M giant. The numbers of common envelope 
systems formed by two-body encounters that lead to total 
envelope expulsion will therefore be dependent on the rem- 
nant fraction in the GC population (see §2). To be sure 
that most common envelope systems will expel the enve- 
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Figure 9. Envelope-ejection efficiency, as a function of impactor 
mass [Mq], for which the envelopes of our giants may be expelled 
before the impactor sinks deep enough into the giant that the 
impactor or the core fills its Roche lobe or the density of the gi- 
ant becomes comparable to that of the impactor. The lower box 
represents the 8 Mq giant and the upper box represents the 2 Mq 
giant. For each giant, the upper-most (solid) line corresponds to 
Roche-lobe filling by ZAMS impactors. The lower, solid line rep- 
resents Roche-lobe filling by stellar remnants; the broken region 
represents remnants that, if they exist, must be produced in a 
binary since their solitary progenitors live longer than a Hubble 
time. The almost vertical drop in the remnant curve corresponds 
to the Chandrasekhar mass limit (~ 1.44 Mq), above which the 
remnants are neutron stars (or black holes for even larger rem- 
nant masses) and, effectively, are point masses requiring very low 
a to eject the envelope. The long-dashed line represents Roche- 
lobe filling of the giant's core and the dot-dashed line represents 
the efficiency at which the ZAMS star falls deep enough into the 
envelope that the local density matches its own. 



lope requires a population with a large abundance of stellar 
remnants or massive main-sequence stars with few low-mass 
main-sequence stars; such a population may exist in the GC 
if the star- formation low- mass cut-off is high (e.g. > 0.8 Mq) 
or if the GC contains an extra population of neutron stars 
and black holes (Morris 1993). 

The violence of this method of forming common enve- 
lope systems may partly assist in lowering the impactor- 
mass limit for full envelope expulsion. Although the en- 
ergy deposited in the envelope to make the impactor bound 
will dissipate on the envelope's thermal time-scale 
of a few-hundred years (perhaps somewhat shorter than 
common envelope evolution), the mass loss, however, is per- 
manent. Both effects tend to lower the envelope's binding 
energy. In both cases, though, the effects are strongest for 
massive impactors (which are likely to eject the envelope 
anyway) and are weakest for low-mass impactors. 



Figure 10. Time-scales [years] for the creation of bound systems 
in the galactic centre. The lower graph represents time-scales for 
the 8 Mq giant and the upper graph corresponds to the 2 Mq 
giant. In each case, the solid line represents 1 Mq impactors, 
the short-dashed line shows time-scales for 2Mq impactors (this 
curve is partly concealed by the solid line for the upper figure), the 
long-dashed line denotes the 4 Mq impactors and the dot-dashed 
line corresponds to 8Mq impactors. In addition, in the 2Mq gi- 
ant portion the dotted curve represents collisions with 0.6 Mq 
impactors and the long-dashed dotted curve shows the bound- 
system creation time-scales with 1.4 Mq impactors. We have used 
the velocity dispersion and density profile of Genzel et al. (1996). 
For both target giants, the two less-bold sets of curves (which 
point to low time-scales at low galactocentric radii) correspond 
to impactor masses of 1 Mq (solid line) and 8 Mq (dot-dashed 
line), but we have replaced the density of Genzel et al. (1996) by 
the p oc r~ 18 profile (see §2.4). 



8 DISCUSSION: TIME-SCALES FOR 
FORMING COMMON ENVELOPE 
SYSTEMS 

The time-scale for a giant to undergo a two-body collision 
that forms a bound system is given by 

1 



ri(r)/ce(r) 



(12) 



where / cc (?") is the fraction of collisions that form bound 
systems and Vi(r) is the collision rate per giant, given by 
Equation If all such bound giant-impactor pairs develop 
into common envelope systems (see §6) which then evolve to 
expel the envelope, r ce is then a lower limit for the destruc- 
tion of our giants in the GC, excluding the rare collisions 
in which the instantaneous mass loss approached the total 
envelope mass (§6). We minimise r co by assuming that, for 
a given impactor mass, M2, all stars in the GC (other than 
our giant) are of mass M2. This produces Figure |l(j where 
we see that for both giants r cc > 10 9 years, even if n oc r -1 ' 8 . 
This time-scale is longer than the lifetimes of stars of masses 
<; 2Mq. Given the rarity of high mass-loss encounters, we 
conclude that the paucity of luminous red-giant stars is un- 
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likely to be due to two-body encounters involving AGB stars 
of mass ~ 2 — 8M0, whether via dramatic, instantaneous 
mass loss or via forming a common envelope system (after 
a collision) which then evolves to expel the envelope. This 
conclusion does not exclude the possibility of collisions in- 
volving pre-giant stars. We also note that collisions involving 
binaries may be another way to destroy giant or pre-giant 
stars (see §9 and Davies et al. 1998). 

The broad minima in Figure nO are due to competition 
between Fi(r) and f co (r). The former is high at low galacto- 
centric radii and decreases outwards, following the trends of 
the velocity dispersion and density profile. f ce (r) is highest 
when the velocity dispersion is lowest (§6.2). The outcome is 
a minimum at about r — 0.3 pc, coincident ally just outside 
the region of the observed paucity. 

Another feature is the narrow spread in r co with im- 
pactor mass. Note that the curves for the 2 Mq giant en- 
countering 2Mq, 4Mq and 8M0 impactors may be un- 
derestimating t cc somewhat [§6.3]). If r cc is independent of 
Ma, then / ce oc M 2 because r cc oc l/(ri/ co ) oc Af 2 /(p/ cc ). 
Including corrections for uncertainties in / ce (f) (see §6.3) 
shows that r cc varies by less than a factor of 2.5 with M2 
across a wide range of radii. Clearly r cc is only weakly de- 
pendent on M2. 



9 FURTHER EXPLANATIONS FOR THE 
OBSERVED DEPLETION OF GIANTS 

9.1 Eccentric orbits 

In the calculations of time-scales for two-body collisions in- 
volving giant stars we have assumed that the stars are on 
circular orbits, i.e. constant galactocentric radii, and so, for 
a given star, the collision rate is uniform in time. A ther- 
malised stellar system, however, has a distribution of ec- 
centricities, and so a number of giant stars will be on orbits 
that take them very deep into the GC. If the stellar number- 
density profile is cusp-like (see §2.4) then the collision rate 
in the deep core will be very high, suggesting that giants on 
such orbits could be destroyed via collisions when near their 
perigees. The time spent by the giants near the perigee is, 
however, short, and preliminary investigation suggests that 
unless the cusp-like density profile (p oc r -1,8 ) assumed here 
significantly underestimates the actual density profile, then 
the number of devastating collisions involving giants will be 
insufficient to account for the observed paucity, even if the 
impactors are massive black holes rather than solar-mass 
stars. 

Main-sequence stars as well as giants can be on eccentric 
orbits that take them deep into the GC, and so possibly the 
progenitors of the giants may be destroyed in collisions. Fig- 
ure pi suggests that collisions involving main-sequence stars 
could become frequent in the deep GC if the number-density 
of stars is cusp-like. Considering main-sequence stars on ec- 
centric orbits of semi-major axes of about 0.12 pc, the num- 
bers of collisions that they will undergo in the deep GC is 
too small to account for the total absence of bright giants as 
observed by Genzel et al. (1996). It is worth noting, however, 
that the longevity of low-mass main-sequence stars indige- 
nous to very low galactocentric radii allows them to undergo 
several collisions. 



9.2 Mass segregation 

Stars of mass greater than the local mean mass will tend to 
sink deeper into the potential well via dynamical friction: 
mass segregation will occur. Mass-segregation occurs on a 
time-scale similar to the relaxation time, which for a star of 
given mass M, is given by: 



T re l = 



0.34cr(r) 



(13) 



G 2 Mp(r)lnA 

(see Binney & Tremaine 1987), and thus depends on the lo- 
cal values of density, velocity dispersion and A, the Coulomb 
potential. As r re i may be a function of radius in the GC, so 
also will be the mass-segregation rate, which may be esti- 
mated as (Hut et al. 1992): 

r —. (14) 

dt T-rel(r) 

Thus the time taken for a star's orbit to segregate from 
radius n to rt is 



0.34 

G 2 M 



a dr 
rplnA 



(15) 



For mass segregation to be of interest to the collisional de- 
struction of stars, we require r sog < r m , where r m is the 
lifetime of a star of mass M. If a given star sinks signifi- 
cantly during its lifetime, then the increase in collisions that 
it will undergo in high-density regions could become impor- 
tant. Disregarding the contribution from A, r scg oc 1/M; 
stars of higher mass sink more quickly than those of lower 
mass. However, r m decreases strongly with Al, so that for 
8M stars sinking from radii of 0.1- 0.2 pc, the increase in 
the numbers of collisions is too small to account for the ab- 
sence of the bright giants (Genzel et al. 1996), if such giants 
were 8Mq stars. Stars of lower masses, say 2Mq and 4Mq, 
live sufficiently long for significant mass segregation to oc- 
cur, as Figure [ll] shows. The cusp-like density profile gives 
shorter segregation times than the profile of Genzel et al. 
(1996) because the slope of the former tends to offset the 
increasing Keplerian (a oc r -0 ' 5 ) velocity dispersion as ra- 
dius decreases, whereas the latter profile is essentially flat at 
low radius and so the increasing velocity dispersion actually 
gives rise to lower sinking rates as radius decreases. 

Although the flatter density profile of Genzel et al. 
(1996) may allow 2Mq and 4Mq stars to sink some dis- 
tance during their lifetimes, the collision rates do not in- 
crease strongly enough to make a difference compared to 
the case of a star that remains at a constant galactocentric 
radius of about 0.12 pc; this is a consequence of the profile's 
flatness (Figure [j). The cusp-like profile, however, not only 
allows stars to sink more deeply into the GC but also pro- 
vides a strongly increasing collision rate. Assuming the cusp- 
like profile, during its lifetime a 4 Mq star may sink down to 
about 0.02 pc if it were created at 0.12pc. By virtue of their 
longevity, stars of mass 2Mq may sink to within 0.01 pc in 
the cusp- like profile, even if they were born at 0.17pc. Whilst 
mass segregating inwards, both stars will undergo the ma- 
jority of collisions whilst on the giant branch, i.e. during the 
latter stages of their lives after they had had time to sink 
some distance. A giant star that has mass segregated from a 
certain galactocentric radius will have more collisions than 
an unsegregated giant star on an eccentric orbit because 
the former will spend all its life at low galactocentric radii 
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radius 



Figure 11. Mass segregation time-scales [in years] for stars as a 
function of galactocentric radius [parsecs] . The time measured is 
the time taken to sink inwards from an initial radius of 0.12 pc 
(heavy lines) and from 0.17pc (light lines) for stars of mass 2Mq 
(continuous lines) and 4Mq (short-dashed lines). In each case, we 
have assumed the usual velocity dispersion profile [Equation (g|)] 
and two values of the density profile: that of Genzel et al. (1996) 
and the cusp-like profile (see §2.4); the former gives rise to a 
longer segregation time-scale and thus, for a given star, is always 
the higher line. The horizontal dotted lines indicate the lifetimes 
of 2Mq stars (upper line) and 4Mq stars (lower line). 

whereas the latter is only present in the deep GC when at 
periastron. A segregated giant on an eccentric orbit will, of 
course, have more collisions than in the other two cases. 

Assuming that all other stars were solar-type, from 
Equation (||) we calculate that a 4Mq star indigenous to 
0.02 pc would undergo ~ 80 collisions, and a 2Mq star on 
a circular orbit at 0.01 pc would undergo ~ 400 collisions. 
Figure [l(] shows that at such low radii, the time-scale for a 
2Mq giant to undergo a collision that leads to a common- 
envelope system is longer than the stellar lifetime; this is 
very likely also true for a 4Mq giant. Thus the formation 
of bound systems cannot destroy enough of the giants to 
account for the paucity. 

The mass lost from the envelope drops with increas- 
ing impact speed, and so collisions at very low galactocen- 
tric radii will involve lower mass losses than the (small) 
quantities reported in §6. Individually, therefore, solar mass 
impactors cannot destroy the giants, but it may be possi- 
ble that the giants could be progressively stripped of their 
envelopes over the many collisions that they would suf- 
fer. Impacts of speeds germane to the very low-radius GC 
(Voo <; 800km/s) with solar-mass impactors will only re- 
move small amounts of envelope mass from either giant, and 
so it is unclear whether solar-mass impactors may be effec- 
tive in destroying giants. Impactors of greater mass may be 
more effective. For a target giant in a sea comprised solely 
of 8Mq impactors (black holes), the number of collisions 
will be reduced by a factor of eight compared to solar-type 



impactors. Cumulative mass loss via several collisions may 
remove most of the giants' envelopes. Alternatively, if a gi- 
ant is likely to undergo several collisions with very massive 
impactors, there may be a significant probability of a very 
low impact-parameter collision. This raises the possibility 
of the core being pulled out of the envelope by the massive 
impactor, thereby destroying the giant. An important prob- 
lem with this high-mass impactor scenario is that 2 Mq stars 
generally cannot sink as far into the GC if the local mean 
mass increases from our assumed 1 Mq to 8 Mq . 

Collisions involving mass-segregating main-sequence 
stars may be important in the destruction of stars before 
they evolve into giants. Unlike the giant phase, the main- 
sequence stars will undergo collisions distributed over the 
range of radii through which the stars sink. For the main- 
sequence counterparts of the two giants considered above, 
the upper limit of the numbers collisions suffered is about 
30 collisions for the 2 Mq star and about two collisions for 
the 4 Mq star, for solar-type impactors and assuming a cusp- 
like density profile. Collisions between main-sequence stars 
at the high speeds involved may lead to a significant amount 
of mass loss and possibly a merger of the two stars, depend- 
ing on the impact parameter and mass ratio of the colliders 
(Benz & Hills 1987, 1992; Lai, Rasio & Shapiro 1993). It 
is quite possible that potential giants may be destroyed via 
mass loss in collisions whilst on the main-sequence, if the 
deep GC density profile is cusp-like. The merging of the two 
main-sequence stars into a higher-mass star has the effect of 
removing two potential lower-mass giants from the galactic 
centre whilst adding a higher-mass star that will eventually 
evolve into a massive giant. If such a process occurs, then 
it is inconsistent with the paucity observed by Genzel et 
al. (1996) because massive giants are extremely luminous 
whereas the observations show an absence of such stars in 
the deep GC. 

Another consideration must be the process of resonant 
dynamical friction (see §9.3), which may take place deep 
in the GC. As the mass-segregation time-scale is similar to 
the relaxation time-scale, resonant dynamical friction may 
give rise to segregation time-scales shorter than the nor- 
mal segregation time by a factor ~ 7M st /Mbh (see Rauch & 
Tremaine 1996), where M st is the enclosed stellar mass inside 
a given radius. This process should not occur for a system in 
which the orbits are isotropic (Rauch & Tremaine 1996). In 
a spherically symmetric, well-mixed distribution of orbits, 
however, the density profile cannot decrease inwards. Obvi- 
ously, the observed paucity of the brightest red giants means 
that the density distribution of those stars does decrease in- 
wardly within the central ~ 0.2 pc, thereby implying that 
the brightest giants' orbits are not well-mixed or isotropic. 
However, resonant dynamical friction allows orbits to drift 
in angular momentum space, rather than in energy space 
and thus their semi-major axes would remain constant. 



9.3 Tidal disruption by a central, massive black 
hole 

The galactic centre likely harbours a 2 x 10 6 Mq black hole 
(Eckart & Genzel 1997). Stars will be disrupted if they ap- 
proach the black hole to within the tidal radius, which is 
given by (Hills 1975); 
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Rt 



( 6Mbh\ 



1/3 



(16) 



where Mbh = 2 x 10 M© and pi is the mean density of 
the star. We find that R t = 7.4 x 1CT 4 pc = 3 x 10 5 R 
for a star of mass 8 Mq and radius 243 R© . For a giant of 
mass 2Mq and radius 94 R©, R t = 4.6 x 10 -4 pc. In other 
words, if these giants had orbital semi-major axes of 0.1 pc, 
they would require eccentricities > 0.993 for them to be 
disrupted. Assuming isotropic orbits, the probability, P t , of 
a star being on an orbit that takes it inside the tidal radius 
is given by (e.g. Sigurdsson & Rees 1997): 

m 

where d is the semi-major axis of the star's orbit. For the 
star of mass 8M and radius 243 R Q , P,~7x 10" 3 . The 
observed paucity of giants may be explained by tidal dis- 
ruption only if, on time-scales shorter than the giants' life- 
times, the orbits of all the giants are scattered into the phase 
space occupied by the loss cone. The time-scale for such scat- 
tering to occur is usually taken to be the relaxation time 
[Equation (fL3|)], but nearly Keplerian potentials (such as 
that near the central black hole) can allow faster resonant 
relaxation to occur (Rauch & Tremaine 1996). This pro- 
cess is important within the galactocentric radius at which 
Monc = M st /(M st + Mbh) = 0.1 (Rauch & Ingalls 1998), 
where M s t is the total enclosed stellar mass. Within the re- 
gion of missing bright giants (<J 0.2 pc), ^i cnc satisfies this 
condition, assuming either the cusp-like density profile or 
that of Genzel et al. (1996) (see §2.4). The resonant relax- 
ation time-scale is given by Rauch & Tremaine (1996): 

Trcs = — Trel- (,lo) 

JWbh 

Well inside the region in which the bright giants are de- 
pleted, say at 0.1 pc, r ros ~ 4 x 10 8 years, assuming a mean 
stellar mass of ~ 1 Mq and assuming the density profile of 
Genzel et al. 1996 (the cusp-like profile actually gives rise to 
longer resonant relaxation time-scales as the enclosed stellar 
mass is higher) . Given the small size of P t and also that the 
resonant relaxation time is longer than the giant phases of 
the stars (see Table ^), we conclude that tidal disruption via 
relaxation cannot account for the paucity. 

Giants may also enter the loss cone via precession if 
the GC stellar distribution (and thus potential) were flat- 
tened from spherical, as the total angular momentum of a 
giant in such a potential is not conserved. If the flattening is 
along, say, the ^-direction, then two of the star's angular mo- 
mentum components (J x and J y ) are not conserved whereas 
J z is conserved. For giant stars on orbits of d — 0.1 pc to 
pass within R t of the black hole, the required values of J z 
represent only a small fraction of the full range available. 
Thus precession into the loss cone in a flattened asymmetric 
potential in the GC cannot account for the missing bright 
giants. 

We conclude that tidal disruption cannot account for 
the paucity of the bright giant stars. 



9.4 Binary stars 

Interactions between binary stars and red giants in the GC 
were considered by Davies et al. (1998). Such interactions 



can explain the observed paucity only if the giants are low 
mass (2M© rather than 8Mq) and also if the stars can be 
hit at earlier stages of their evolution. A high binary fraction 
in the GC population is also required. If these conditions are 
satisfied then hard binaries will be more destructive towards 
giants than single stars. For stellar collisions to occur in a 
binary-single star encounter, one expects i? m in Z~> d. The 
collision rate for such encounters is less dependent on the 
individual size of the interacting stars than are two-body en- 
counter rates. Thus main-sequence stars, which are smaller 
but longer lived than giants, will play a greater role in 2+1 
body encounters than in single-single collisions. If a sufficient 
number of collisions occur in 2+1 body encounters then the 
precursors of the bright-giant stars may be destroyed before 
they ascend the giant branch. This will be investigated in 
future work. 

For sufficiently tight binaries, the size of the primary's 
Roche lobe limits the primary's size, or may even lead to 
the premature end of the primary's evolution. The binary 
separation, dhs, for which single cluster members have just 
enough kinetic energy to break up the binary is given by: 



dhs 



GM a M b (M a + M b + Mc 
M c (M a + M b ) 



V 2 



(19) 



(e.g. Davies 1995) where Voo is the relative velocity between 
the binary and the single star; M a , Mb and M c are respec- 
tively the masses of the two binary components and the mass 
of a third-body intruder. Deep in the GC, the mass of the 
central black hole dominates, forcing stars to follow Kep- 
lerian orbits (see Eckart & Genzel 1997). As the velocities 
and velocity dispersions of the stars are independent of the 
their masses, the relative velocity dispersion between a bi- 
nary ab and single star c is therefore <r ab c = <r ab + o 2 = 2<r 2 
and the mean-square velocity is (V^,) = 3<r ab c = 6a 2 . Sub- 
stituting (V^j) for Vx in Equation (|l9|), we find for o(r) 
appropriate for the GC, d^s is indeed small (<J 2Rq for 
M a = Mb = M c — IMq). The time-scale for two objects 
to pass within, say, 100 R© of each other deep in the GC 
is ~ 10 7 years. With the exception of the tightest binaries, 
most binaries will have been ionised before the primary could 
evolve off the main sequence. We therefore expect most of 
the GC giant stars to be single and so the observed paucity 
of bright giants cannot be explained in terms of giants' evo- 
lution being restricted in tight binaries. Another problem 
with this scenario is that it does not explain why the bright 
giants disappear rather precipitously at about 0.2 pc. 



10 CONCLUSIONS 

Our chief conclusions may be summarised as follows: 

(i) The mass losses that result from single-single stellar 
encounters involving red giants in the galactic centre are 
almost always low (^ 0.1 Mq). Only a very small number 
of collisions will involve sufficient mass loss to destroy the 
giant. Such collisions require low f? m i n (^ 30 R©) and pref- 
erentially high-mass impactors (<; 8M©). Such high mass 
loss collisions are thus unlikely to account for the observed 
paucity of bright red giants. 

(ii) Mass loss was found to be an increasing function of 
impactor mass and decreasing functions of i? m in and Voo 
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reflecting that gravitational forces rather than shock mech- 
anisms dominated the collisions. 

(iii) Tidal capture of impactors by giant stars is unlikely 
to occur in the galactic centre. Most collisions involved 
the impactor passing through the giants and remaining un- 
bound. This is due to the high velocity dispersion associated 
with the galactic centre. 

(iv) The fraction of the impactor's orbital energy (AE/E) 
deposited in the giant increased with impactor mass and de- 
creased with increasing i? m in and Voo- For sufficiently small 
-Rmin and Voo the impactor became bound to the giant, im- 
mediately forming a common envelope system or an eccen- 
tric binary system which would rapidly decay into a common 
envelope system. The fraction of collisions resulting in these 
bound systems was ^ 0.5 at a few pc, decreasing to ~ 1CP 5 
at very low radii, depending upon the impactor mass. At the 
radius at which the luminous giants are missing (~ 0.2 pc), 
the fraction was ~ 10~ 2 - 10~ 3 . 

(v) We showed that many of these common envelope sys- 
tems formed would lead to total envelope ejection for the 
2 Mq giant and also for the 8 Mq giant if the efficiency was 
high or the impactor was a high-mass main-sequence star 
or a stellar remnant. However, even if all common envelope 
systems destroy the giant, the time-scale for collisions which 
result in these bound systems was at least the entire lifetimes 
of stars of masses <; 2 M Q and so could not account for the 
paucity of the giants. 

(vi) The time-scale for creating bound systems was only 
weakly dependent of the impactor mass and so can reveal 
little about the nature of the undetectable stellar population 
of the galactic centre. 

(vii) Giants (or their precursors) within the GC may be 
destroyed in collisions involving high mass losses if their or- 
bits are sufficiently eccentric to take the stars deeper into 
the galactic nucleus where the collision rate is higher. This 
requires, however, that we currently significantly underesti- 
mate the stellar density deep in the GC. 

(viii) Dynamical friction and resonant dynamical friction 
may allow stars, originally on wide orbits at about 0.2 pc 
(the outer radius at which bright giants are observed to be 
missing), to travel deep into the galactic nucleus, thereby 
becoming subject to a greatly increased collision rate. This 
study does not rule out the possibility of the bright giants 
(or their precursors) being destroyed in such collisions. This 
process is only significant, however, if the galactic-centre has 
a cusp- like density profile, rather than the observed flattened 
profile. 

(ix) The absence of the giant stars is not due to tidal dis- 
ruption around a central black hole of mass 2.5 x 10 6 Mq, 
nor does it seem likely to be due to single-single stellar col- 
lisions involving giants. To explain the paucity in terms of 
stellar collisions requires investigation into the collisional ef- 
fects of binary stars upon pre-giant stars. Hard binaries in 
the galactic centre may well have separations too small for 
the primary to evolve into a giant. This also, however, means 
that most binaries will have been ionised and so most giants 
will be single stars (unless, perhaps, they exchange into a 
binary). It also does not seem to explain why the numbers 
of bright giants plummet rapidly at ~ 0.2 pc. 
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